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The present article is about tiie statistical mechanics of non-trivial field configurations. The non- 
trivial fields arise from the negative sign of the commutators and the anticommutators of the bosonic 
and fermionic field excitations. These kinds of fields were previously studied by Pauli and Lee and 
Wick. The thermal distribution function of the above mentioned fields are calculated in the article 
and using the thermal distribution functions the energy density, pressure and entropy density of the 
non-trivial field configurations are found out. The results match exactly with a previous calculation 
done by Fornal et. al. for higher derivative Lee-Wick theories showing a deeper similarity with 
the earlier work. It is assumed that such kinds of non-trivial fields may have existed in the early 
{3J[), universe and may have some cosmological relevance. 
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I. INTRODUCTION 



pi i' For the last couple of decades there has been a series of papers related to the Lee- Wick kind of field 

theories. The Lee- Wick construction [l|, |2] originally was motivated to tackle the problem of infinities in 
flj^' quantum field theories. By including massive partners of the photon and other fermions present in the theory 

r^ I of Quantum electrodynamics Lee and Wick could produce a theory virhich was devoid of divergences. Later 

there has been a flurry of activity concerning the connection of higher derivative theories and the Lee- Wick 
construction 0,|j|- ^H these theories assumed the existence of some partners of the Standard model particles. 
The partners of the Standard model fields arose out of suitable redefinitions of fields in a higher derivative 
r^ ' version of the Standard model Lagrangian. In these cases the propagators of the partners of the Standard 

^f) ', model fields behaved more like the propagators one would get if they had a Lee- Wick like theory. There has 

lO ■ been at least one attempt Q to use the concepts of these Lee- Wick constructions in cosmology where the 

authors were able to show the bouncing nature of the universe whose energy is dominated by the energies of 
[~^ ■ a scalar field and its Lee- Wick partner. In Ref. [6[ the authors tried to formulate a possible thermodynamic 

f^ ' theory of particles which includes the Lee- Wick partners using a method of statistical field theory previously 

^D . formulated by Dashcn, Ma and Bernstein in Ref. Q. 

In 1943, Pauli wrote a paper on an unusual field theory [8|, where the bosons appearing in the theory were 
accompanied by partners, whose excitation modes were quantized with the wrong sign of the commutators. 
. J, ■ The idea of proposing such an unusual theory came from an urge to suppress ultraviolet divergences appearing 

rS I in normal quantum field theories. Most of the work of Lee and Wick [l|, Q utilize the concepts of Pauli, 

jrt ' and his field excitations which will be called Pauli oscillators in this paper. Pauli's paper did not have any 

ingredients from a higher derivative theory. 

There has not been any attempt to produce a statistical field theory based on Pauli's work up till now. 
One of the important outcomes of this article is to show that the results obtained on the thermodynamics of 
a Lee- Wick ghost infested universe, as done in [6|, matches exactly with the thermodynamic quantities of an 
universe full of Pauli oscillators. We have included the fermionic sector also, where the basic quantization 
prescription was given by Lee and Wick [2]. In this regard it should be noted that unlike the fermions in 
Lee and Wick's theory the fermion quantization presented in this article is simpler and the fermions do not 
have a complex mass. These facts may have some bearing on the unitarity of the theory but presently the 
theory predicts thermodynamic results which arc in agreement with the results obtained in Ref. [6| . 

The previous attempt to theorize the thermodynamic properties of a universe which is constituted by 
the Standard model particles and their possible Lee- Wick partners gave interesting results which may have 
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very important consequences in cosmology. In the present work we use the basic commutation relations 
and number operator language as used in the original paper of Lee and Wick, and Pauli [g, and derive 
the expressions of the thermal distribution functions of the Lee- Wick ghosts or the Pauli oscillators which 
accompany the normal particles present in the Standard model. The distribution functions of such highly 
non-trivial field configurations are presented for the first time in this article. The thermal distribution 
function of the unusual field configurations are different in nature from those of the standard Bose-Einstein 
or Fermi-Dirac distributions. Using the new distribution functions we can reproduce all the thermodynamic 
quantities previously calculated by Fornal, Grinstcin and Wise in Ref. [6|. 

The material in the present article is presented in the following manner. The next section discusses about 
the technique to find out the thermal distribution function of the Lee- Wick partners. In section IIIII the 
energy density, pressure and entropy density of a gas comprising of elementary particles and their unusual 
field partners are calculated using the thermal distribution functions. The last section IIVI summarizes the 
important points of the article. 



II. DISTRIBUTION FUNCTIONS OF THE PAULI OSCILLATOR EXCITATIONS AND 

LEE- WICK FIELDS 

A free quantum field theory can be thought of to be made up of an infinite number of linear harmonic 
oscillators in the momentum space oscillating with frequencies e(p) = \/pM-rn^ where p is the 3-momentum 
of the oscillator excitations which can take infinite values, and m is the mass of the excitations over the 
vacuum. For each of the oscillators, corresponding to any p, there corresponds an annihilation operator a(p) 
and a creation operator a(p). The operator a(p) is defined as a(p) = rj^^a^ {p)ri , where t] is the metric on 
the Hilbert space of the oscillators which is -1-1 for most of the cases, but need not be so in general. The 
oscillators can be bosonic, if the only non- vanishing commutators are of the form [a(p) , a(k)] = (5^(p — k) , 
or fermionic if the only non- vanishing anti-commutators arc of the form {a(p) , a(k)} = (5'^(p — k) . In the 
above cases the metric on the Hilbert spaces of the bosonic or fermionic oscillators are both -1-1. But the 
above mentioned commutation relations do not exhaust all the possibilities of quantizing the oscillators. 
Instead of the usual bosonic commutation, if we assume 



[a(p),a(k)] = -,53(p-k), 
then also we get oscillator like spectrum where the eigenvalues of the number operator [2] : 

iV(p) = -a(p)a(p) 



(1) 



(2) 



are all positive, and the minimum value of A^(p) is zero. But in this case the metric rj on the Hilbert space of 
the oscillators become indefinite, which implies it is not always -1-1 but its general form is 77 = (— 1)~^(p^. In 
the present case we can define a vacuum satisfying A''(p)|0) ~ 0. For such field excitations with momentum 
p the Hamiltonian of a single oscillator can be written as 



^(P) = 2<P) HpHp) + a(p)a(p)] = -<P) 



1 



^(p) + o^ (0) 



(3) 



The additive term — ^ e(p) (5'^(0) is the zero point energy which can be ignored in further calculations. In the 
present case the Hamiltonian of the above field excitations are completely equivalent to the negative energy 
Hamiltonian of the Pauli oscillator in Ref. Q . Lee and Wick in Ref. [l|, 0] did show that it is indeed possible 
to produce a sensible unitary quantum field theory with such kind of excitations. In this article we are more 
interested about the thermal distribution of the Lee- Wick or the Pauli excitations. Some one interested in 
the formal structure of such kind of theories can refer to the original work in Ref. [l|, [^ . Due to the presence 
of on-shell excitations the thermal vacuum becomes \Q) = |n(pi), n(p2), • • • } where n(pi) is the number of 
excitations carrying momentum pi. The action of the number operator on such a vacuum is [2| 



iV(p)|fl) = n(p)|f]), 



(4) 



where n{p) is the number of particles with momentum p present in the thermal vacuum. In this analysis 
we will consider non-interacting real scalar field for which the chemical potential ^ = 0. In the present case 
with a Hamiltonian of the form as given in Eq. ([3]) the single particle partition function will be, 
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n(p)=0 



„/3n(p)e(p) 



(5) 



where /3 



From the last equation it is seen if we are finding the single particle partition function for the 



Lee- Wick partner of a Standard model boson then the series representing 2;^ does not converge for /3 > 0. 
Consequently we regularize the last expression by cutting off the summation for a finite value of n(p) as: 
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where Af is a phenomenological cut-off which can be made indefinitely big at the end of the calculation. 

Next we calculate the thermal distribution function of the field excitations, which behave like Pauli oscil- 
lator excitations, from the expression of the single cell partition function of the field excitations as given in 
Eq. ^. In conventional statistical mechanics we can find the single cell distribution function via 
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(7) 



V,/3 



where /i is an auxiliary chemical potential whose exact nature is not important for our purpose. In presence 
of an auxiliary chemical potential the single particle partition function can be written as 
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Applying conventional methods, the distribution function can also be written as 

/b(p) 
which comes out to be, 
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Now setting the auxiliary chemical potential to be zero, y^t = 0, we get the distribution function of the fields 
whose excitation spectra is given by the Pauli oscillators as 
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This is the distribution function of the fields whose excitations are described by creation and annihilation 
operators which satisfy Eq. ([Ij. These are not the fields which appear in the Standard model of particle 
physics. The distribution function as plotted in Fig. [T] shows that the average excitation per energy level is 
negative definite. Obviously these systems describe a physical theory which is non-trivial and the negative 
sign of the distribution is only meaningful when compared with the positive definite distribution of the 
normal Standard model bosons. In general these kind of distributions will produce negative energy density 
and pressure but once these energy density and pressure is added with the positive energy density and 
pressure of the Standard model bosons we get a net positive energy density and pressure. The important 
point to notice about the distributions is that there is no pile up of quantas near e(p) = as is in the case 
of the Bosc-Einstcin distribution. The reason being that the Pauli excitation spectra distribution has two 
infinite spikes as the energy tends to zero and they cancel each other near the origin. 

If we assume the anticommutators of the creation and the annihilation operators which define the fermionic 
excitations as 



a'(p)=«s(p) = 0, {a,(p),a,.(k)} 



"5' (P - k) , 



(12) 



where s, s' may be some internal quantum numbers, then we can proceed in a similar way as done before and 
calculate the thermal distribution of these excitations. If we stick to the definition of the number operator 
as given in Eq. ([2]), then in this case the possible eigenvalues of the number operator are simply and 1. 
The Hamiltonian of a single oscillator is 



-ff(p) = 2^*^^^ [as(p)as(p) - as{v)as{v)] = -e(p) 



N{p)-\6\0) 



(13) 
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FIG. 1: The plot of the distribution function as given in Eq. fTT}. The topmost curve is for a normal Bose-Einstein 
distribution and the lower two curves correspond for /b(p) for the ultraviolet cutoff M — 100 and 200. The inverse 
temperature in all the cases is .OlGeV"^ and the energy e(p) is in GeV. 

where 2 e(p) '^''(0) is the zero point energy. In the above equation we have assumed that the Hamiltonian 
is independent of s. If we use this oscillator Hamiltonian to calculate the single particle partition function, 
neglecting the zero point energies, there is no problem regarding the convergence of the series. The single 
particle partition function for the anticommuting fields turns out to be 



n(p)=0 



;/3n(p){«(p)-M} = 1 + g/3{<^(p)-M} 



(14) 



where /x is an auxiliary chemical potential. Now applying the formula in Eq. ^ and setting fi = a,t the 
end we get the distribution function for the Lee- Wick ghost partners or the Pauli excitations of the Standard 
model fermions as: 



/f(p) = — 



1 



,flt(p) + 1 



(15) 



Unlike the previous case, in the present scenario the distribution function has no dependence on the dimen- 
sionless regulator M. 

Here we point out that the fermionic degrees of freedom are not exactly the same as presented in the 
paper by Lee and Wick. There they assumed that the fermionic ghosts mix with each other and more over 
they have complex masses. The description of the fermionic sector in this article is more in line with the 
fermionic extension of Pauli's work. It is to be noted that although the Standard model particles arc proper 
bosons or fermions whose thermal distributions are given by the Bose-Einstein or Fermi-Dirac distributions, 
in the present theory the bosons (fields obeying commutation relations) or the fermions (fields obeying 
anti-commutation relations) do not follow the Standard Bose-Einstein or Fermi-Dirac statistics. 

^From the nature of the distribution functions, for both the bosonic and fermionic excitations, it is seen 
that the average number of particles in an energy cell e(p) is negative. This is a very counter intuitive 
result. But once we observe both the bosonic and fermionic Hamiltonians we notice that the energy of any 
single excitation of these fields are negative. The negative energy is not a result of any particular choice of 
potential or interaction of the fields but rather an artifact of the way these fields are quantized. The vacuum 
defined is not stable and there exists much less energetic states than the vacuum itself. These kinds of fields 
are unstable. The maximum energy of the field configurations can be zero. As the distribution functions are 
interpreted as average number of excitations of the field quantas in any energy level e(p), which is positive, 
we immediately notice a contradiction. In the present case the field configurations cannot afford to have 
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FIG. 2: The plot of the distribution function as given in Eq. (|15|1 . The topmost curve is for a normal Fermi-Dirac 
distribution at /3 = .OlGeV"^ and the lower two curves correspond for /f(p) for /? values .OlGeV^^ and .lGeV~^ 
and the energy e(p) is in GeV. 

quantum excitations in positive energy states. Consequently the only way there can be some positive energy 
out of these systems is by removing some of the oscillator modes, or particles, from these systems. As a 
result the negative distribution functions indicate that the field configurations must be de-excited to get to 
positive energy states i.e. the Lee- Wick partners decay into the Standard model particles. 



III. ENERGY DENSITY, PRESSURE AND ENTROPY DENSITY FROM THE 

DISTRIBUTION FUNCTION. 



A. The bosonic case 
1 

To calculate the relevant thermodynamic quantities for the bosonic field from a statistical mechanical point 
of view we will employ Eq. (|lip . The energy density can be calculated using the following known equation 



//•OO 
e(p)fB{e)d'p^^J^ <p)fB{e)\p\'d\p\, 



(16) 
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where g stands for any intrinsic degree of freedom of the particle. For a rclativistic excitation e^ = p'^ + m^ 
where m is the mass of the bosonic excitations. Changing the integration variable from |p| to e one gets 
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In the above integral it is assumed that | p | ^ to and to have a closed integral the lower limit of the integral 
is assumed to be zero. In the extreme rclativistic limit the system temperature T ^ to. Both of the integrals 
can only be done when /3 < 0, and in that case the result of the last integral is 



g 

"27r2 



^4r4 TO^TT^T^ 



15 



12 



_g_ fn^_rrM^ 
2^ V15Af3 12A/ 



(18) 



Analytically continuing the above result for /? > and taking Af — > oo we sec that for normal temperatures 
the energy density for extreme rclativistic excitations of the bosonic fields is of the following form 
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(19) 



As expected, the energy density turns out to be negative for the excitations in this case. The pressure density 
of the bosonic field excitations can be found out from 



P-^I ^/bW^^P = ^ f ^/B(s)d|p| . (20) 
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FoUowing similar steps as in the case of the energy density, it is seen that the pressure of extremely relavistic 
excitations of the bosonic fields turns out to be 

The entropy density of the bosonic fields is simply given by 

These values of the energy density, pressure and entropy density exactly match the corresponding values 
calculated for the Lee- Wick partners, in [6[, in a different way. In [6| the authors were trying to formulate 
thermodynamics for a higher-derivative theory. The higher derivative theory was converted into standard 
theory (theory up to a second derivative) with the introduction of Lee- Wick partners. The authors in the 
previous work did not quantize the system explicitly but were working with the form of the propagators of 
the Lee- Wick partners. 

If we assume that in the early universe for each bosonic degrees of freedom in Standard model there exist 
a corresponding bosonic degree of freedom whose creation and annihilation operators are quantized as in 
Eq. ([IJ then the net energy density, pressure and entropy density of the early universe turns out to be 

gm^T'^ grn^T'^ gm^T , ^ 

PB = PSM + P= — , PB = PSM +P= — , SB = SSM + S= , (23) 

which arc all positive as expected. Here the energy density, pressure density and entropy density for the 
Standard model bosonic particles are psM = g ^gT , PsM = 9 ^ on ^^^^ ^sm = g ^'^.g respectively (9|. 

B. The fermionic case 

In this subsection we apply Eq. (|15p to find the energy density, pressure and entropy density of the 
fermionic excitations. In this case the distribution function do not have any dependence on the regulator M . 
For relativistic excitations the integrals which give the energy density and pressure for the fermionic case 
are exactly similar with the bosonic case except that now we have to use the distribution for the fermions. 
The integrals can be easily done, granted /3 < 0, but the results can be analytically continued for positive 
temperatures. The results in this case are listed below. The energy density, pressure and entropy density of 
the Lee- Wick partners are as follows: 

P^-.r^^-^V (25) 

^ -^ V 720 48 y ' ^ ' 

' = ' [-Im - ^ j ■ (26) 

The energy density and pressure quoted above was derived in [y] in a different way for the special case 
of g = 2. If we assume that to each unusual fermionic degree of freedom there corresponds one standard 
fermionic degree from the Standard model, then the total fermionic contribution is 

gm'^T'^ gm'^T'^ gm'^T , . 

PF = PSM + P= ^g , PF = PSM +P= ^g , SF == SSM + S = ^^ (27) 

which are also all positive. Here the energy density, pressure density and entropy density for the Standard 

T 2rT74: 'T-_.2 rri4: 7-— 2 rfi3 ^^T 

Model fermionic particles are psM = g 240 ' Psm = 9 720 ^'^^ ■^sm = ffn^cp respectively [9[. 

It is worth pointing out here that higher derivative theories of fermions require two auxiliary Lee- Wick 
partners (one left-handed and the other right-handed) to eliminate the higher derivative terms. In that case 
the Lee- Wick degrees of freedom exceeds the one of its Standard model partner yielding negative energy, 
pressure and entropy density [lO|. This issue is yet to be resolved. 



IV. DISCUSSION AND CONCLUSION 

In this article we have studied a system of bosonic and fermionic fields, whose excitations modes are 
quantized with the negative sign of the commutator or anticommutator. The bosonic part of the theory was 
well studied before by Pauli Q and the fermionic part is studied previously by Lee and Wick [2]. These 
unusual quantization process naturally produces field configurations whose total energy is negative. The 
negative energy of the field configuration is not due to any particular form of the potential but solely an 
outcome of the quantization process. The vacuum of the theory is not the state with the lowest energy, it is 
rather the state with the maximal energy making the field configuration unstable. The bosonic and fermionic 
degrees of freedom do still follow commutation and anticommutation relations and specifically the fermionic 
fields still follow the Pauli exclusion principle. In this article the emphasize had been on the calculation 
of energy density, pressure and entropy density of the unusual filed configurations. To calculate the above 
mentioned thermodynamic quantities one requires to have a statistical mechanics of the field excitations. 
One encounters the difficulty of a diverging sum when calculating the single particle partition function of 
the bosonic fields. Keeping to conventional ways, where the temperature of the system is positive definite, 
the partition function can only be summed when one uses a ultraviolet cutoff. The distribution function 
calculated from the partition function turns out to be negative definite, which is a nontrivial result. The 
negative nature of the distribution function implies that there must be an average loss of particles in any 
positive energy level. This fact can be understood by noticing that the system can only have positive energy 
by loosing negative energy quanta and the negative sign of the distribution function implies such a condition. 

It is shown in the article that the entropy density of the non-trivial field configurations turn out to 
be negative. This is a serious issue as by its very nature entropy is always positive. The problem with the 
negative entropy of the Lee- Wick partners is avoided if we include the entropy of the Standard model particles 
also. The total bosonic or fermionic entropy turns out to be positive. We presume that the negative entropy 
of the Pauli excitations and Lee- Wick partners show that they cannot exist alone, the thermodynamic system 
is only complete if the Standard model particles are also included. 

The energy density, pressure calculated from the distribution functions of the unusual fields discussed in 
this article are exactly the same as calculated by Fornal, Grinstein and Wise in Q . The derivation of the 
new distribution functions and using them to show the apparent connection between the Pauli oscillator 
thermodynamics and the thermodynamics of the Lee- Wick partners are the main motivations for this work. 
Apparently the work of the previous authors and the work in the present article has a common conceptual 
similarity. The higher derivative theories and their lower derivative Lee- Wick partner infested theories were 
motivated to kill the divergences in a quantum field theory. The original work of Pauli [8| also aimed 
to produce a quantum field theory devoid of divergences. We strongly believe that the work of Pauli and 
Grinstein, O'Connell and Wise have an underlying equivalence although the two theories arises from different 
perspectives. In short, the present article acts as a bridge between higher derivative theories and normal 
theories with deformed quantization conditions. We hope that similar to the Lee- Wick theory, there can 
be a theory of the early universe where all the normal fields are accompanied by Pauli bosons or fermions 
which are quantized with a negative commutator or anticommutator bracket. This kind of a theory has the 
potential of damping the infinite divergences which crop up in normal quantum field theories. The unusual 
field configurations discussed in this article may have been present in the very early universe and, like the 
Lee- Wick ghosts, have decayed to other particles as time went on. 
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this article. 
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